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1 Introduction 

In this paper, we consider the following nonlinear heat equation 



u t --§z\li (x, t) u x ] + f{u) = fi(x, t), < x < 1, < t < T, 



(1.1) 



associated with conditions 



u x (Q,t) = h u(0,t) + g (t) 



u x (l,t) = hxu(l,t) +gi(t), 



(1.2) 



and initial condition 



u(x,0) = Uq(x), 



(1.3) 



where uq, p,, /, /i, go, Pi are given functions satisfying conditions, which will be specified 
later, and ho, hi > are given constants, with ho + h\ > 0. 

The conditions fll.2)l are commonly known as Dirichlet - Robin conditions. They 
connect Dirichlet and Neumann conditions. Theses conditions arise from the effect of 
excess inert electrolytes in an electrochemical system through perturbation analysis 
(0, 0, 0, 0). 

The governing equations (II. 2p are the equation usually used in a diffusion, con- 
vection, migration transport system with electrochemical reactions occurring at the 
boundary electrodes and submitted to non linear constraints. 

In electrochemistry, the oxidation-reduction reactions producing the current is mod- 
eled by a non linear elliptic boundary value problem, linearization of which gives the 
Dirichlet - Robin conditions (0). Theses conditions also appear in the response of an 
electrochemical thin film, such as separation in a micro - battery. His analyze is made 
by solving the Poisson - Nernst - Planck equation subject to boundary conditions 
appropriate (Dirichlet - Robin conditions) for an electrolytic cell (0). 

The paper consists of six sections. In Section 2, we present some preliminaries. 
Using the Faedo - Galerkin method and the compactness method, in Section 3, we 
establish the existence of a unique weak solution of the problem (11.11) - (11. 3p on (0, T), 
for every T > 0. In section 4, we prove that if the initial condition is bounded, then so 
is the solution. In section 5, we study asymptotic behavior of the solution as t — » +oo. 
In section 6 we give numerical results. 

2 Preliminaries 

Put fl = (0, 1), Qt = fl x (0,T). We will omit the definitions of the usual 
function spaces and denote them by the notations L p = 1/(0,), H m = H m (ft) . Let 
(■,■) be either the scalar product in L 2 or the dual pairing of a continuous linear 
functional and an element of a function space. The notation 1 1 ■ 1 1 stands for the norm 
in L 2 and we denote by || ■ \ \x the norm in the Banach space X. We call X' the dual 
space of X. We denote L p (0,T; X), 1 < p < oo the Banach space of real functions 
u : (0,T) — > X measurable, such that ||w||lp(o,T;X) < +oo, with 

(/ T |Kt)||^) 1/P , if l<P<oo, 
ess sup||ii(t) 1 1 x> if p = oo. 

0<t<T 

Let u(t), u'(t) = u t (t) = u(t), u x (t) = \/u(t), u xx (t) = Au(t), denote u(x,t), 
§ (x,t), g(:M), fit(M), respectively. 

On H 1 we shall use the following norms \\v\\ Hl = (\\v\\ 2 + Hi^H 2 ) 1 ^ , IMI, = 

(^) + IHI 2 ) 1/2 ,* = o,i. 



\ u \ \LP(0,T;X) 



2 



Let /i € C° (Qt) , with fi(x,t) > /i > 0, for all (x,t) 6 and the constants 
h , hi > 0, with Hq + hi > 0, we consider a familly of symmetric bilinear forms 
{a(t; -)}o<t<T on H l x H 1 as follows 

a(t; u, v) = f fi(x, t)u x (x)v x (x)dx + h [i (0, t) u(Q)v(0) + h\{i (1, t) u(l)v(l) 

= (fi{t)u x , v x ) + h fi (0, t) u{0)v{0) + hin (1, t) u(l)v(l), for all u, v e H 1 , < t < T. 

(2.1) 

Then we have the following lemmas. 

Lemma 2.1. The imbedding H 1 C°([0, 1]) is compact and 



\v 



c°(n) 



< v2 ||f ||^i , for all v E H 1 , 

< V2 \\v\l , for all v e H\i = 0, 1. 



(2.2) 



Lemma 2.2. Let /i G C° , wift yu(x,t) > /xq > 0, /or a// £ <5^5 anc ^ ^ e 

constants h , hi > 0, with h + hi > 0. Then, the symmetric bilinear form a(t\ -, ■) is 
continuous on H 1 x i/ 1 and coercive on H , i.e., 



(2.3) 



(i) |a(t; m, f)| < a T Hvll^i , 

(m) a(t; > a , 

/or a// m,d6 if 1 , < t < T, where ay = (1 + 2/i + 2/ii) sup //(x, £), and 

(x,t)eQ T 

{/i min{/i , §}, /to > 0, /ii > 0, 
(2.4) 
/i min{/ii, |}, /ii > 0, h > 0. 

The proofs of these lemmas are straightforward. We shall omit the details. 
Remark 2.1. It follows from (12 .2[) that on if 1 , v i — > \\v\\ Hl and v i — >■ ||f 1^ are 
two equivalent norms satisfying 

75 < \\ v \\i < v 7 ^ Hull^i , for all v e H\ i = 0, 1. (2.5) 
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3 The existence and uniqueness theorem 

We make the following assumptions: 
(Hi) h >0 and hi > 0, with h + h x > 0, 
(H 2 ) u eL 2 , 
(H 3 ) g , gi eW^(0,T), 

(H 4 ) ^^([0,1] x [0,T]), fi{x,t) >//o>0, V(z,t) G [0,1] x [0,T], 
(ff 5 ) heL\0,T;L 2 ), 

(H 6 ) f G C°(R) satisfies the condition, there exist positive constants Ci, C{, C2 and p > 1, 
(i) >C x \u\ p -C' x , 
(2i) |/(«)| < C 2 (l + K" 1 ), for all !i£l. 

The weak formulation of the initial boundary valued (II. ip - (jl.3p can then be given 
in the following manner: Find defined in the open set (0, T) such that u(t) satisfies 
the following variational problem 

£(u(t),v) + a(t, u(t),v) + (f(u),v) = (h(t),v) - fi (0, t) g o (t)v(0) - p (1, t) </i(*Ml), 

(3.1) 

Vf G if 1 , and the initial condition 

u(0) = M - (3.2) 
We then have the following theorem. 

Theorem 3.1. Let T > and (Hi) — (H e ) hold. Then, there exists a weak solution 
u of problem (II. ip - (11.31) such that 

u G L 2 (0,T; H 1 ) nL°°(0,T; L 2 ), 

(3.3) 

tu G L°°(0, T; tf 1 ), iut G L 2 (0, T; L 2 ). 

Furthermore, if f satisfies the following condition, in addition, 

(H 7 ) (y - z) (f(y) - f(z)) >-5\y- z\ 2 , for all y, z G M, with 5 > 0, 

then the solution is unique. 

Proof. The proof consists of several steps. 

Step 1: The Faedo - Galerkin approximation (introduced by Lions [5]). 
Let {wj} be a denumerable base of H 1 . We find the approximate solution of the 
problem (II. ip - (11.31) in the form 

Um(t) = J2T=l C mj(t)Wj, (3.4) 
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where the coefficients c m j satisfy the system of linear differential equations 
(u' m (t),Wj} + a(t; u m (t), Wj) + (f(u m (t)), Wj) 

= (fi(t), Wj ) - fi(0,t) g (t) Wj (0) - fi(l,t) gx^Wjil), 1 < j <m, 

, M m (0) = U 0m , 

(3.5) 

where 

uom = J2f=i a mjWj u strongly in L 2 . (3.6) 

It is clear that for each m there exists a solution u m (t) in form (13.41) which satisfies 
(13.51) and (I3.6P almost everywhere on < t < T m for some T m , < T m < T. The 
following estimates allow one to take T m = T for all m. 

Step 2. A priori estimates. 

a) The first estimate. Multiplying the j th equation of (13. 5p by c m j (t) and summing 
up with respect to j, afterwards, integrating by parts with respect to the time variable 
from to t, we get after some rearrangements 

||M m (t)|| 2 + 2 f*a(s;u m (s),u m (s))ds + 2 f Q (f(u m (s)),u m (s))ds 

= ||w 0m || 2 + 2 fo(h{s),u m (s))ds (3.7) 

-2 f fj l (0,s)g (s)u m (0,s)ds-2 f* (1, s) p 1 (s)« fn (l, s)ds. 

By u 0m — > u strongly in L 2 , we have 

IKm|| 2 < Co, for all m, (3.8) 

where Co always indicates a bound depending on m . 

By the assumptions (H e , (i)), and using the inequalities (I2.2p . (12.31) . and with (3 > 0, 
we estimate without difficulty the following terms in (13. 7p as follows 

2 f a(s;u m (s),u m (s))ds > 2a f*\\u m (s)\\ 2 H ids, (3.9) 

2 H (f(u m (s)), u m (s))ds > 2d J* \\u m (s)\\ p LP ds - 2TC[, (3.10) 
2/o(/i(s),u m (s))ds < ||/i|| £ i (0) t;l2) + /o \\u m (s)\\ 2 ds, (3.11) 

-2 Jj fj, (0, s) g (s)u m (0, s)ds < 2^2 ||//|| Loo(Qt) \\g \\ Loo Jj ||w m (s) ds 

< | T IHIl°°(Q t ) Ibollioo +Pfo \\Um{s)f H x ds, 

(3-12) 

-2 Jj /i (1, s) #i(s)w m (l, s)ds < 2^2 ||/i|| ioo(Qr) \\ gi || Loo J * ||u m (s) || H1 ds 

< f T IMlL°(Q T ) IMIl- +/?/o ||Wm(s)|lii ^, 

(3.13) 
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for all > 0. Hence, it follows from ([37]) - (13TT3]) that 

||w m (t)|| 2 + 2(a - /?) /o ||t*m(s)||^i ds + 2Ci /„* ||« m (s)||^ P 

< Co + + ||/i|| Ll(0iT;i2) + /J ||/i(s)|| \\u m (s)\\ 2 ds (3.14) 

+ f IMIi,~(Q T ) (|M|L» + IblHioc) • 

Choosing (3 = |a , we deduce from ( I3.14p . that 

5 m (t) < 4 1} + /o 4 2) ( s )^m(s)^, (3-15) 

where 

S m (t) = \\u m (t)\\ + a J Q || ds + 2d £ || (s) ||^ p (is, 

4 1} = Co + 2TC1 + ||M| Ll(0iT;L2) + ±T \\n\\ 2 Laa{QT) (|bo||^ + IMlLc) , (3.16) 

I c?oo = ilAOOII. 4 2) e^(0,T). 

By the Gronwall's lemma, we obtain from (I3.15p . that 

S m (t) < 4 X) exp (/ '4 2) (^s) < C T , (3.17) 

for all m £ N, for all i, < t < T m < T, i.e., T m = T, where Ct always indicates a 
bound depending on T. 

b) The second estimate. Multiplying the j th equation of the system (13. 5p by t 2 c' m j(t) 
and summing up with respect to j, we have 

\\tu' m (t) f + t 2 a(t- u m (t), u' m (t)) + (tf(u m (t)), tu'Jt)} 

= (th{t)M' m {t)) - 1 2 » (0, t) g (t)u' m (0, t) - t 2 ix (1, t) gi{t)u' m (l, t). 

(3.18) 

First, we need the following lemmas. 
Lemma 3.2. 

(i) § (t; u, v) = (// (■, t) u x , v x ) + ho/j,' (0, £) u(0)v(0) + h lf i' (1, t) w(l)«(l), /or a// u,v E H\ 



(ii) || (t; m, t>) | < ar ||w||#i \\v\\ h i , /or a// w, v £ if 1 , 

(t),M m (t)), 

(3.19) 



at' 

(hi) |a(t; u m (t),u m (t)) = 2a(t; u m (t) , u' m (t)) + (t; u m (t) , u m (t)) , 



where a T = (1 + 2/i + 2^) sup //' (x, t) . 

(x-,t)e[0,l]x[0,T] 



Lemma 3.3. Put A = , m = J^ o |/(|/)| d?/, and /(z) = f* f(y)dy, 



Z £ 



Then we have 



-mo < f(z) < C 2 (\z\ + i |^D, e R. (3.20) 

The proofs of these lemmas are straightforward. We shall omit the details. ■ 
By (I3.19IK . we rewrite (I3.18P as follows 

2 \\tu' m (t)\\ 2 + f t a(t;tu m (t),tu m (t)) + 2(tf(u m (t)),tu' m (t)) 

= 2ta(t;u m (t),u m (t)) + § (t;tu m (t),tu m (t)) +2(tf 1 (t),tu' m (t)) ( 3 - 21 ) 
-21* n (0, t) g o (t)u' m (0, t) - 2t 2 n (1, t) gi(t)u' m (l, t). 
Integrating (I3.2ip . we get 
2 fa \\ su ' m ( s )\\ 2 ds + a(t;tu m (t),tu m (t)) + 2 f*(sf(u m (s)),su' m (s)}ds 



2 f* sa(s; u m (s),u m (s))ds + J dt (s; su m (s), su m 



We shall estimate the terms of (I3.22p as follows. 

a(t;tu m (t),tu m (t)) > a \\tu m (t)\\ 2 Hl , 

2 Jo(sf(u m (s)), su' m (s))ds = 2 jl s 2 ds£ £ dx j^ s) f(y)dy 

= 2 f* s 2 dsf s f(um(x, s))dx 



(3.22) 
(3.23) 



2 Jo fs ( s2 Jo f( u m(x } s))dxj - 2s fl f(u m (x, s))dx 
2t2 Jo f( u m(x, t))dx - 4 f* sds Jg 1 f(um(x, s))dx 



ds 



> -2T 2 m -4C 2 f s \\u m {s)\\ L i + \ \\u m {s 



\LP 



> -2T 2 m - ATC 2 



^ \\ U m\\L°°(0,T;L 2 ) + p 2Ci^ m ^ 



2j Q sa(s;u m (s),u m (s))ds < 2Ta T J \\ U m {S) \\ H i ds < 2Ta T ±S m {t) < Q 
f Q ff (s; su m (s), su m (s))ds < a T f Q \\su m (s)\\ 2 H1 ds < T 2 a T /„ 



■T, 

u m {s)\\w ds 



< T 2 a T ±S m (t) < C T , 



ds 

> -C 7 
(3.24) 

(3.25) 
(3.26) 
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2 fo( S fl( S )i SU 'm( S )) ds < 2 Jo \\ S h( S )\\ \\ SU 'm( S )\\ ds < Jo \\ S h( S )f ds + Jo \\ SU 'm( S )\\ 2 ds 



(3.27) 



<r 2 /o T |l/l(«)l| 2 ^ + /oll™m( S )H 2 ^ 

<C T + f*\\su' m (s)\\ 2 ds. 
By using integration by parts, it follows that 
-2 So S V (°> s ) 9o(s)u' m (0, s)ds 

-2t 2 y (0, t) g (t)u m {0, t) + 2 jj [s 2 n (0, s) g (s)}' u m {0, s)ds 

< 2V2t 2 |MU (Qt) \\g \\ L ~ \\u m (t)\\ Hl + 2V2f*\[s 2 v(0,s)g (s)}'\ \\u m (s)\\ Hl ds 

< INI£~ Wt) IMIi« + /? ||*«m(*)fc + 2v/2 /J | [s 2 /i (0, s) g (s)]'\ \\u m (s)\\ m ds 
<^C T + (3 \\tu m {t)f Hl + 2V2 f* I [s 2 /i (0, s) g {s)}'\ \\u m (s)\\ H1 ds. 

On the other hand 

I [s 2 /i (0, s) g (s)]'\ = \2sfx (0, s) g {s) + s 2 [// (0, s) g (s) + y. (0, s) g' {s)} 



(3.28) 



< 2s 



=(Q T ) \\9o\\ L oo + s 2 M\ C \Q T ) [\\9o\\ L <« + lfo(' 



(3.29) 



< a M CH Q T ) [(2 + T) Ibollxoc + T < sCWo(s), 

where 

Cr=|Hlci(5 r) [(2 + T)|| flb || i0O +T], ^00 = 1 + 1^001, VoeL^T). (3.30) 
Hence, we deduce from ( I3.28p . (13.291) . that 

-2j Q t s 2 f i(0,s)go(s)u' m (0,s)ds < \C T + /3 \\tu m (t)\\ 2 m + 2\[2Ct So 4>o(s) \\su m (s)\\ H1 ds 

<±C T + (3 \\tu m (t)\\ 2 Hl + 2C 2 SoMs) d * + So \\su m (s)\\ 2 m ds 



< (1 + \)C T + /3 II^WH^ + £ Vo(s) 



for all /3 > 0. 
Similarly 

-2 £ s 2 /i (1, s) gi{s)u' m {l, s)ds < (1 + i)C T + \\tu m (t)\\ 2 Hl + f* Vx(s) \\su m (s 



(3.31) 



(3.32) 



for all j3 > 0, where 

C , T=||^|| C1( Q T) [(2 + T)||( 7l || ioo +T], ^(s) = l + |^(s)|, Vie^(0,T). (3.33) 
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It follows from - dS2ZD, f lCT]) and that 

/ * ||sM' m (s)|| 2 cis + a ||tu m (t)||^i 

< (6 + f)C T + 2/3 11^(011^ + Jo Ms) \\su m (s)f Hl ds (3.34) 
+ SoMs) \\su m (s)f Hl ds. 
Choosing 2/3 = ^a , we deduce from (13.341) . that 

X m (t) < + SoCT\s)X m (s)ds, (3.35) 

where 

' X m {t) = \\tu m {t)\\l 1 + Sl\\su' m {s)\\ 2 d Sl 

< C? = (l + ±)(G + ^)C T , (3.36) 

k c?(s) = (1 + i) (Vo(«) + , 4 2) e ^(o,r). 

By the Gronwall's lemma, we obtain from (13. 35 p . that 

+ Jo ll*<(*)l| 2 ^ < ^?exp (/ T cf (s)ete) < C T , (3.37) 

for all m £ N, for all t £ [0, T], VT > 0, where Ct always indicates a bound depending 
on T. 

Step 3. The limiting process. 

By (13.161) . (I3.17P and (13.371) we deduce that, there exists a subsequence of {u m }, 
still denoted by {u m } such that 

u m — y u in L°°(0,T;L 2 ) weak*, 

u m — > u in L 2 (0, T; if 1 ) weak, 

< tu m -±tu in L oc (0,T;.H' 1 ) weak*, (3.38) 

(tu m )' — > (t-u)' in L 2 (Q T ) weak, 

Mm — ^ u in L p (Qt) weak. 

Using a compactness lemma ([5], Lions, p. 57) applied to (j3.38p ^/i. we can extract 
from the sequence {w m } a subsequence still denotes by {u m }, such that 

tw m — > tu strongly in L 2 {Qt). (3.39) 

By the Riesz- Fischer theorem, we can extract from {u m } a subsequence still de- 
noted by {u m }, such that 

u m (x, t) ->■ u(x, t) a.e. (x,t) in Q T = (0, 1) x (0, T). (3.40) 
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Because / is continuous, then 

f(u m (x,t)) ->• f(u(x,t)) a.e. (x,t) in Q T = (0, 1) x (0, T). (3.41) 

On the other hand, by (H e , ii), it follows from (I3.16p . (I3.17P that 

\\f(u m )\\ LP , {QT) <C T , (3.42) 

where Cy is a constant independent of m. 

We shall now require the following lemma, the proof of which can be found in [5] . 
Lemma 3.4. Let Q be a bounded open set of 1R and G m , G G L q (Q), 1 < q < oo, 



\\G r 



and 



< C, where C is a constant independent of m, 



G m ^rG a.e. (x,t) in Q. 



Then 



G m G in L q (Q) weakly M 

Applying Lemma 3.4 with N = 2, q = p' , G m = f(u m ), G = f(u), we deduce from 
(EOT]) . 0332D, that 

f{u m ) -> /(m) in L p '{Q T ) weakly. (3.43) 
Passing to the limit in (13. 5 p by (13. 6p . (I3.38p . (I3.43p . we have satisfying the equation 



dt 



[u(t),v) + a(t, u(t),v) + (f(u), v) 



(fi(t),v) - fi (0, t) g (t)v(Q) - /x (1, t) gi (t)v(l), Vv e F 1 , (3.44) 



k u(0) = u . 



Step 4- Uniqueness of the solutions. 
First, we shall need the following Lemma. 

Lemma 3.5. Let u be the weak solution of the following problem 



{ u t - £[fi{x,t)u x } = f(x,t), 0<x<l, 0<t<T, 
u x (0, t) - h o u(0, t) = u x (l, t) + h lU (l, t) = 0, 
u(x, 0) = 0, 

u G L 2 (0, T; H 1 ) n L°°(0, T; L 2 ) n L P (Q T ), 
{ tu G L°°(0, T; H 1 ), tu t eL 2 (Q T ). 



(3.45) 



Then 



(3.46) 



IK*)l| 2 + 2 /o a(s,ix(s),ii(s))ds = 2 f*(f(s),u(s))ds. 

The lemma 3.5 is a slight improvement of a lemma used in [T] ( see also Lions's 
book 0).B 
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Now, we will prove the uniqueness of the solutions. Assume now that (H 7 ) is 
satisfied. 

Let Mi and u 2 be two weak solutions of (11. ip - (11. 3p . Then u = u\ — u 2 is a weak 
solution of the following problem (I3.45P with the right hand side function replaced by 
f(x,t) = —f(ui) + f{u 2 ). Using Lemma 3.5 we have equality 

\\u(t) f + 2 jj a(s, u(s),u(s))ds = -2 - f(u 2 ),u(s))ds. (3.47) 

Using the monotonicity of /(?/) + we obtain 

/„'(/(«!) - /(«„), u(a))<fa > -5fi \\u(s)\\ 2 ds. (3.48) 

It follows from (E3ZD, ([33SD that 

\\u(t)\\ 2 + 2a jt\\u(s)f Hl ds<25j*\\u(s)\\ 2 ds. (3.49) 

By the Gronwall's Lemma that u — 0. 
Therefore, Theorem 3.1 is proved. ■ 



4 The boundedness of the solution 

We now turn to the boundness of the solutions. For this purpose, we shall 
make of the following assumptions 

(H[) h >0 and hi > 0, 

(H' 2 ) u GL°°, 

(H' 5 ) fi G L 2 (Q T ), fx(x,t)<0, a.e. (x,t) e Q T , 

(H'q) f G C°(R) satisfies the assumptions (Hq) , (H?) , and 

uf{u) > 0, Wu G K, |«| > H^oIIloo ■ 
We then have the following theorem. 

Theorem 4.1. Let (H[), (H 2 ), (H 3 ), (H 4 ), (H' 5 ), (H' 6 ) hold. Then the unique weak 
solution of the initial and boundary value problem (11.11) - (11.31) . as given by theorem 
3.1, belongs to L°°(Q T ). 

Furthermore, we have also 

\\ u \\l^(Qt) — max { II^oIIloo , ^ IIs'qII^^t) ' ~ti Ibiili^^.T)} (^-i) 

Remark 4.1. Assumption (-f/^) is both physically and mathematically natural in 
the study of partial differential equation of the kind of (II. ip - (II. 3p . by means of the 
maximum principle. 

Proof of Theorem 4.1. First, let us assume that 

u (x) < M, a.e., x G Q, and max |i Hfl-oLc^T) , ^ IMIz^ct)} < M. (4.2) 
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Then z = u — M satisfies the initial and boundary value 
f z t -£[^(x,t)z x } + f(z + M) =h(x,t), 0<x <1, 0<t<T, 

z x (0,t) = ho[z(p,t)+M\+9o(t), -z x (l,t) = h 1 [z(l,t)+M\+g 1 (t), (4.3) 
z(x, 0) = Uq(x) — M. 

Multiplying equation (I4.3p 1 by v, for v G H 1 integrating by parts with respect 
to variable x and taking into account boundary condition (I4.3)) q. one has after some 
rearrangements 

Jq 1 z t vdx + J fi (x, t) z x v x dx + fi (0, t) [h (z(0, t) + M) + g (t) } v(0) 

+fi(l,t) [h 1 (z(l,t) + M)+g 1 (t)]v(l) 

+ J 1 /(z + M)vdx = J fi(x, i)vdx } for all v E H 1 . 

(4.4) 

Noticing from assumption we deduce that the solution of the initial and bound- 
ary value problem JTTJ - jT3) belongs to L 2 (0,T; /7 1 ) nL°°(0, T; L 2 )nL p (Q T ), so that 
we are allowed to take v = z + = ^(\z\ + z) in (14. 4p . Thus, it follows that 

2 t 2; + rfx + /„* /x (x, *) z^+dz + /j (0, t) [/i (z(°> t) + M) + g (t) ] z + (0, t) 

+li (1, *) t) + M) + 9l (t) } z+(l, t) (4.5) 

+ Jq 1 /(z + M)z + dx = Jq 1 /i(x, t)z + dx. 

Hence 

|| ||^(t)|| 2 + a(t,^(t),^(t)) + / 1 /(^ + + M)^rfx = $lh{x,t)z+dx 

-ii (0, t) (/ioM + g (t)) z+(0, t)-fi (1, t) {hM + 9l (t)) z+(l,t) < 0. 

(4.6) 

since 

M > max{^ \\g \\ Loo , 7^ ML-} and 

J z t z + dx = f z>0 (z + ) t z + dx = ~ Jo z>0 |^ + | ^=2dt/ol z+ l ^ X = 2~dt • 

(4.7) 

and on the domain z > we have z + = z and z x = (z + ) x . 

On the other hand, by the assumption (H' 2 ) and the inequality (12. 3p . we obtain 

a(t,z+(t),z+(t))>a \\z+(t)\\ 2 m . (4.8) 

Using the monotonicity of f(z) + 5z and (H 7 ) we obtain 

/o /(^ + + M)z+dx = J a [f(z+ + M) - /(M)] z+cfe + J,, 1 f{M)z + dx 

> -6 J* \z + \ 2 dx + Jo f(M)z + dx > -6f* \z + \ 2 dx = -S\\z^ 



(4.9) 
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Hence, it follows from (J£S)L (Q]| . that 



2 tit 



^(t)r + a |k + (t)||^<5||2+(t)|r (4.10) 

Integrating (I4.10p . we get 

||^(t)|| 2 <||^(0)|| 2 + 25/ c ;||^( s )|| 2 rf S . (4.11) 

Since z + (0) = (u(x,0)—M) + = (u (x) — M) + = 0, hence, using Gronwall's Lemma, 
we obtain ||z + (t)|| 2 = 0. Thus z + = and u(x, t) < M, for a.e. (x,t) G Qt- 

The case —M < u (x), a.e., iffl, and M > max H^oII^^t) , 7^ llfi'i lli°°(o,T) } 
can be dealt with, in the same manner as above, by considering z = u + M and 
z~ = h(\z\ — z), we also obtain z~ = and hence u(x,t) > —M, for a.e. (x,t) G Qt- 

From all above, one obtains \u(x,t)\ < M, a.e. (x,t) G Qt, i.e., 

IMIz- (0t ) < M, (4.12) 

for all M > max |||m || loo , £ ||5o|| x <-(o,r) > 7TT IMIz°°(o,t)} ' 
This implies ( 14. ip . Theorem 4.1 is proved. ■ 

5 Asymptotic behavior of the solution as t +00. 

In this part, let T > 0, (Hi) — (Hi) hold. Then, there exists a unique solution u of 
problem (II. ip - (11.31) such that 

u G L 2 (0,T; H 1 ) fl L°°(0, T; L 2 ) D L P (Q T ), 

tue L^iO^H 1 ), tv! G L 2 (Q T ). 

We shall study asymptotic behavior of the solution u(t) as t — > +00. 
We make the following supplementary assumptions on the functions /1 (x, t) , f x (x, t) , 
9i(t), 92(t). 



(m) g ,gieW 



1,1/ 



fi G ^([0, 1] x R+), MM) > ^0 > 0, V(x,t) G [0, 1] x R+, 

(ifg) There exist the positive constants Ci, 71, gooo, gioo and the functions 
G C^fO, 1]), /ioo G L 2 , such that 
(«) |Po(*) -Pood <Cie-^, Vt>0, 
(ii) - fl'iool <Cie-^, Vt>0, 

(m) ||yu (t) - AioolLoo < Cie"^ 1 *, Vt > 0, //^(x) > fi > 0, Vx G [0, 1], 

||A(t) -Aooll <c ie -^, vt >o. 
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First, we consider the following stationary problem 

d_ 

dx 



■§- [/ioo (x) U x \ + f(u) = floo(x), < X < 1, 



(5.1) 

u x {0) = h u(0) + g 0ixi , - u x (l) = hiu(l) + gf loo . 

The weak solution of problem (15.11) is obtained from the following variational prob- 
lem. 

Find £ H 1 such that 

^OO (^OO J ^ ) + (f{Uoo),v) = (fioo,v) ~ ^(0)^000^(0) - Hoo(l)giooV(l), (5.2) 

for all v G H 1 , where 

a^iu.v) = f Q fi 00 (x)u x (x)v x (x)dx + /i ^oo(0)«(0)w(0) + /i 1 // 0O (l)w(l)w(l) 

= (HooUx, v x ) + h fioo(fyu(0)v(0) + /ii^ 0O (l)u(l)v(l), for all u,v E H 1 . 

(5.3) 

We then have the following theorem. 

Theorem 5.1. Let (Hq) , (H^) — (Hq) hold. Then there exists a solution Moo of 
the variational problem (15.21) such that G H 1 . 

Furthermore, if f satisfies the following condition, in addition, 

{H") f(u) + 5u is nondecreasing with respect to variable u, with < 5 < do- 

Then the solution is unique. 

Proof. Denote by {wj}, j = 1,2, ...an orthonormal basis in the separable Hilbert 
space H l . Put 

m 

Urn = Y.dmjWj, (5.4) 

i=i 

where d m j satisfy the following nonlinear equation system: 

aooiVm, Wj) + (f(y m ),Wj) = (/ioo, Wj) - /ioo (0)5000^(0) ~ l^oo ( l)PlooWj ( 1) , 1 < < TXl. 

(5.5) 

By the Brouwer's lemma (see Lions [5], Lemma 4.3, p. 53), it follows from the 
hypotheses (H e ) , {Hi[) — (Hq) that system (15.41) . (15.51) has a solution y m . 



Multiplying the j equation of system (I5.5P by d m j, then summing up with respect 



■th 

to j, we have 

^oo {ym,y m ) + (f(Vm),Vm) = (fl oo 5 ym) /^oo 

(0)g 0oo y m (0) - ^00(1)01 oo2/m 

(1). (5.6) 

By using the inequality (12.31) and by the hypotheses (i^) ; (#3 ) ~ (^O > we obtain 

+ C1II < Ci + [H/looH + V^d/ioo (0)<?0oo I + |/ioo(l)(?loo|)] II 

(5.7) 

Hence, we deduce from (15 .7p that 

Iz/mllfp ^ C, 

(5.8) 

I ym \\lp ^ C, 
14 



C is a constant independent of m. 

By means of (I5.8P and Lemma 2.1, the sequence {y m } has a subsequence still 
denoted by {y m } such that 

y m -> Moo in H 1 weakly, 

y m —> in L 2 strongly and a.e. in Q, (5.9) 
y m — > Moo in L p weakly. 

On the other hand, by (I5.9IW and (-He) > we have 

/(l/m) -> /(woo) a.e. in ft. (5.10) 
We also deduce from the hypothesis (if 6 ) and from fl5.8j) o that 

Jo l/(U*))r' < 2 p '- 1 Cf'[l + £ \y m {x)\ p dx] < C, (5.11) 

where C is a constant independent of m. 

Applying Lemma 3.4 with N = 1, q = p' , G m = f(y m ), G = /(woo), we deduce 
from f lBTTOjl . (ETTA that 

fiVm) -> /(woo) in L p ' weakly. (5.12) 

Passing to the limit in Eq. (I5.5p . we find without difficulty from (I5.9p . (15.121) that 
Uoo satisfies the equation 

OooKjo, Wj) + (f{Uoo),Wj) = (/loo, Wj) - fioo(0)g 0oo Wj(0) - (l)#looWj (1) . (5.13) 

Equation (15 . 13[) holds for every j = 1, 2, i.e., (15.21) holds. 

The solution of the problem (15.21) is unique; that can be showed using the same 
arguments as in the proof of Theorem 3.1.1 

Now we consider asymptotic behavior of the solution u(t) as t -> +oo. 
We then have the following theorem. 

Theorem 5.2. Let (H t ) , (H 2 ) , (H 6 ) , (H%) - (Hg) , (H!}) hold. Then we have 

\\u(t) - Uoo \\ 2 < (\\u - Uoo f + ip ^ y ) e~ 2 ^, V* > 0, (5.14) 

where 

< 7 < min{7x, a — 5 — 4e}, < As < a Q — 5, 

C > is a a constant independing of t. 

Proof. Put Z m (t) = u m (t) — y m . Let us subtract f)3.5P i with (I5.5P to obtain 

(^IW.Wj) + a(t;u m (t),Wj) - a 00 (y m ,i«j) + (f(u m (t)) - f(y 

m)i w j) 

(fi(t) - /loo, Wj) - [A* (0, t) g Q (t) - /ioo(0)# oo ] Wj(0) 

(5.15) 

- [/i (l,£)#i (t) - fioo{l)g loo }wj(l), 1 < j < m, 



Jm ym- 
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By multiplying (15. 15j) i by c m j(t) — d m j and summing up in j, we obtain 



Id. II 7 
2 dt W^™ 



+ a(t; Z m (t), Z m (t)) + a(t; Vmi Z m (t)) Oqo [y m , Z m (t) ) 
+(/(« m (t))-/(y m ),Z m (t)} 

= - /loo, Z m (t)) - \M (0, t) O (f) " /ioo(0)(?0oc ] Z m (0, t) 

- \fi(l,t)gi(t) -//oo(l)^loo ]Z m (l,t). 



(5.16) 



By the assumptions (i^') — (Hq) , (H7) , and using the inequalities (12.2 j) . (12. 3j) . and 
with e > 0, we estimate without difficulty the following terms in (I5.16P as follows 



a(t; Z m (t), Z m (t)) > do \\Z ri 



(f(u m (t))-f(y m ),Z m (t)) > -5\\ Z u 



> -5\\Z n 



(5.17) 
(5.18) 



(*)) = Mt) P'oo) Vmxi Z mx 

(*)> 

+h (ja (0, t) - ^oo(O)) y m (0)Z m (0, t) (5.19) 

(/i (1, t) - /ioo(l)) y m (l)Z m (l, t); 
Note that ||2/ TO ||#i < C, we obtain from (15. 19j) that 

\a(t; WVmxW \\Z mx {t)\\ 

+2h \\/j,(t) -/ioollioo Ibmllffi \\Z m (t)\\ Hl 

+2hi II// (t) - Aioo|| LO o \\y m \\m \\ Z m(t)\\ H i 

< (1 + 2h + 2/n) C ie -^*C \\Z m {t)\\m < z \\Z m {t)\\ 2 H i + \Ce~ 2 ^- 

\(fl{t) ~ floc,Z m {t))\ < ||/i(t) - / loo || ||Z m (*)|| 



(5.20) 



<c ie -^|l^ 



< £ ||Z m (t)||^ 1 + iCe- 2 ^; 



(5.21) 



\fi (0,t)g (t) - (0)5-000 ]Z m (0,t) 

= - [(/x (0, t) - ^(0)) ffo (f) + ^(0) (g (t) - g 0oo ) } Z m (0, t) 

< \/2||Z m (t)|| ffl \\(i(t) -/ioo|| L oo ||»o|Il«»(R + ) +Moo(0) |^o (*) -^oc 



< v^2||Z m (t)|| H1 |bo|| 



Il°°(k+) 



+ /ioo(0) 



Cie~ 71 * <e||Z n 



^ + ±Ce" 2 ^. 



Similarly 



- [12 (l,t) 9l (t) - ^(Ikioo ] Z m (l,*) ||^m(i) \\m + iCe- 2 ^. 



(5.22) 



(5.23) 
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It follows from flgjgp - (j53Sj) . ( lo^Uj) - and Q, that 

4 ll^(t)|| 2 + 2 (oo - 5 - 4e) \\Z m (t)\\m < fCe" 2 ^. 



c2t 



(5.24) 



Choose £ > and 7 > such that ao — 5 — As > and 7 < min{7i, ao — 5 — 4e}, 
then we have from (15. 24p that 

|2 , n.. 11 ^ ^Ml2 < 8^ e -2 7l ^ 



d* II ^ m 



(5.25) 



-27* 



(5.26) 



+ 2y\\Z, 

Hence, we obtain from (I5.25P that 

\\Z m (t)f< [\\Z m (0)\\ 2 + dfbfi 
Letting m — > +00 in (I5.26j) we obtain 

\u(t) - Uoof < hminf \\u m (t) - y m \\ 2 < (j\u - u^f + ^hfi ) e ~ 27 *, for all < > °. 

(5.27) 

This completes the proof of Theorem 5.2.B 



6 Numerical results 

First, we present some results of numerical comparison of the approximated 
representation of the solution of a nonlinear problem of the type (II. ip - (11. 3p and the 
corresponding exact solution of this problem. 
Let the problem 



where 



u t - u xx + /(it) = fi(x,t), < x < 1, t > 0, 
u x {0, t) = 2u(0, t) + £,,(*), - t) = u(l, t) + ^i(t), 
u(x, 0) = ^0(^)7 

' /i(x,t) = -e x (l + 2e-*) + (l + e-*) p - 1 e( p - 1 ) a; , 
/( u ) = |m| p ~ 2 m, p= |, 
(/„(*) = -1 - e"*, <?i(t) = -2e(l + e-% 



(6.1) 



(6.2) 



k ito (x) = 2e x . 

The exact solution of the problem (16. ip . (16. 2p is u(x, t) = (1 + e~ t )e x . 
To solve numerically the problem (16. ip . (I6.2p . we consider the nonlinear differential 
system for the unknowns = u(x k ,t), Xk = kh, h = 1/N. 

f ^{t) = ±U k -i - j^u k + j^u k+ i - f(u k ) + fi{x k ,t), 
M o = j^K ( M i ~ hg {t)) ,u N = T ^ K (m w _i - hgi{t)) , 
k u k (0) =u (x k ), k = 1,2, N — 1. 
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or 



^t(t) = ^Mfe-i - + ^u k+1 - f(u k ) + £), fc = 2,N-2, 



du 



N- 



u k (0) = u (x k ), k = l,N-l. 

(6.3) 

To solve the nonlinear differential (16. 3p at the time t, we use the following linear 
recursive scheme generated by the nonlinear term f(u k ) : 



= W (St) «?° + ^ - /("M) - it) + h{x u t), 



(n) 



* H*) = - & (t§t) - " mtU)^(^) + 

{ 4 n) (0) =M (x fe ), fc = l,iV-l. 

(6.4) 

The linear differential system ( 16. 4 p is solved by searching the associated eigenvalues 
and eigenfunctions. With a spatial step h — | on the interval [0, 1] and for t € [0,3], 
we have drawn the corresponding approximate surface solution (x,t) — > u(x,t) in 
figure 1, obtained by successive re-initializations in t with a time step At = For 
comparison in figure 2, we have also drawn the exact surface solution (x, t) — > u(x, t). 

Note that, the approximate solution u(x,t) decreases exponentially to Uqq{^x) as 
t tends to infinity, being the unique solution of the corresponding steady state 
problem 

I 3 

— u xx + \u\ 2 u = —e x + e2 x , < x < 1, 

(6.5) 

u x (0) = 2u(0) - 1, -u x (l) = u(l) - 2e. 



18 



approximated solution 




Figure 1. Approximated solution 

exact solution 
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